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We use a deformed QRPA formalism to describe simultaneously the energy distributions of the 
single beta Gamow- Teller strength and the two-neutrino double beta decay matrix elements. Calcu- 
lations are performed in a series of double beta decay partners with A — 48, 76, 82, 96, 100, 116, 128, 
130, 136 and 150, using deformed Woods-Saxon potentials and deformed Skyrme Hartree-Fock mean 
fields. The formalism includes a quasiparticle deformed basis and residual spin-isospin forces in the 
particle-hole and particle-particle channels. We discuss the sensitivity of the parent and daughter 
Gamow- Teller strength distributions in single beta decay, as well as the sensitivity of the double 
beta decay matrix elements to the deformed mean field and to the residual interactions. Nuclear 
deformation is found to be a mechanism of suppression of the two-neutrino double beta decay. The 
double beta decay matrix elements are found to have maximum values for about equal deformations 
of parent and daughter nuclei. They decrease rapidly when differences in deformations increase. We 
remark the importance of a proper simultaneous description of both double beta decay and single 
Gamow- Teller strength distributions. Finally, we conclude that for further progress in the field it 
would be useful to improve and complete the experimental information on the studied Gamow- Teller 
strengths and nuclear deformations. 

PACS: 21.60.Jz, 23.40.Hc, 23.40. Bw 



I. INTRODUCTION 

The recent experimental confirmation of neutrino oscillations has reinforced the interest of nuclear processes in- 
volving neutrinos, see [1] and references therein. Nuclear double beta decay is a rare second order weak interaction 
process that takes place when the transition to the intermediate nucleus is energetically forbidden or highly retarded. 
Two main decay modes are expected in this process. The two neutrino mode, involving the emission of two electrons 
and two neutrinos, and the neutrinoless mode with no neutrino leaving the nucleus. While the first type of process is 
perfectly compatible with the Standard Model, the second one violates lepton number conservation and its observation 
is linked to the existence of a massive Majorana neutrino. For this reason, considerable experimental and theoretical 
effort is being devoted to the study of this process [2] . 

From the theoretical point of view, one particular source of uncertainty is the evaluation of the nuclear matrix 
elements involved in the process. They have to be calculated as accurately as possible to obtain reliable estimates for 
the limits of the double beta decay half-lives. Since the nuclear wave functions and the underlying theory for treating 
the neutrinoless and the two-neutrino modes are similar, the usual procedure is to test first the nuclear structure 
component of the two-neutrino mode against the available experimental information on half-lives. 

The proton-neutron quasiparticle random phase approximation (pnQRPA or QRPA in short) is one of the most 
reliable and widely used microscopic approximations for calculating the correlated wave functions involved in f3 and 
double beta decay processes. The method was first studied in Rcf. [3] to describe the f3 strength and was also 
successfully applied to the description of double beta decay [4] after the inclusion of a particle-particle (pp) residual 
interaction, in addition to the particle-hole (ph) usual channel. Many more extensions of the QRPA method have 
been proposed in the literature, see Ref. [5] and references therein. 

An extension of the pnQRPA method to deal with deformed nuclei was done in Ref. [6] , where a Nilsson potential was 
used to generate single particle orbitals. Subsequent extensions including Woods-Saxon type potentials [7], residual 
interactions in the particle-particle channel [8], selfconsistcnt deformed Hartree-Fock mean fields with consistent 
residual interactions [9] and selfconsistcnt approaches in spherical neutron-rich nuclei [10], can also be found in the 
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literature. In Refs. [9,11], /?-dccay properties were studied on the basis of a deformed selfconsistcnt HF+BCS+QRPA 
calculation with density dependent effective interactions of Skyrme type. A deformed QRPA approach based on a 
phenomenological deformed Woods-Saxon potential was used to calculate the Gamow- Teller strength distributions 
for the two decay branches in double beta decay of ^^Ge [11,12]. 

Recently, the issue of nuclear deformation, which has usually been ignored in the QRPA-like treatments of nearly 
spherical nuclei, was raised in Ref. [12,13]. In Ref. [12] it was found that differences in deformation between the 
initial and final nuclei can have large effects on the double beta decay half-lives. Within the deformed QRPA using 
axially-symmetric symmetric Woods-Saxon single particle basis the particular case of the two-neutrino double beta 
decay (2i^/3/3-decay) of ^^Ge was analyzed [12]. The effect of deformation on the other double beta decay processes of 
experimental interest has not been sufficiently studied yet, see [12,14,15] and references therein. 

In this work we extend these deformed calculations [9,12] by studying first the single (3 branches that build up the 
double beta process and then the 2i^/3/3-decay process itself. We focus on the f3~ Gamow- Teller (GT) transitions of 
the double beta emitters as well as on the /J"*" Gamow- Teller transitions of the daughter nuclei ending up to the same 
intermediate virtual nucleus. The cases considered are those where the two- neutrino double beta decay half- lives have 
been measured, namely 



In sect. II, we present a brief summary containing the basic points in our theoretical description. Section III 
contains the results obtained for the ground state properties of the nuclei mentioned above. In sect. IV we present 
our results for the GT strength distributions and discuss their dependence on the deformed mean field and residual 
interactions. Sect. V contains the results for the two-neutrino double beta decay calculations. The conclusions are 
given in sect. VI. 



In this Section we summarize briefly the theoretical formalism used to describe the Gamow- Teller transitions. More 
details can be found in Refs. [9,12,11]. 

The single particle energies and wave functions are generated from two different methods to construct the deformed 
mean field, which is assumed to be axially symmetric. In one case we start from a deformed Woods-Saxon potential. 
The parameters of this potential are taken from Ref. [16]. The isospin dependence of this parametrization allows one 
to extend it to any mass region. Previous QRPA calculations have shown that it provides a realistic description of 
the ground state properties of deformed nuclei as well as good results on Ml excitations [17] for nuclei in various 
mass regions. The quadrupole deformation of the WS potential is determined by fitting the microscopically calculated 
quadrupole moment to the corresponding experimental value. The hexadecapole deformation is expected to be small 
for these nuclei and we assume it is equal to zero. 

In the other method we follow a selfconsistcnt Hartree-Fock procedure to generate microscopically the deformed 
mean field. This is done with density dependent effective interactions of Skyrme type. Contrary to the previous case, 
the equilibrium deformation of the nucleus is obtained now selfconsistently as the shape that minimizes the energy 
of the nucleus. In this work we present the results obtained with the most common of the Skyrme forces, Sk3 [18] 
although sometimes we also show for comparison results obtained with the force SG2 [19]. 

In both schemes, WS and HF, the single-particle wave functions are expanded in terms of the eigenstates of an 
axially symmetric harmonic oscillator in cylindrical coordinates, using eleven major shells in the expansion. Pairing 
correlations between like nucleons are included similarly in both cases in the BCS approximation with fixed gap 
parameters for protons and neutrons. The gap parameters are determined phenomenologically from the odd-even 
mass differences through a symmetric five term formula involving the experimental binding energies. The values 
obtained from this procedure for the nuclei under consideration can be seen in Table 1. 

The deformed quasiparticle mean field is now complemented with a spin-isospin separable residual interaction that 
contains two parts, an attractive particle-hole and a repulsive particle-particle. The coupling strengths of these forces 
Xc^T ^^"^ ^GT defined as positive. The strength of the ph force is usually determined by adjusting the calculated 
positions of the GT giant resonances to experiment. The strength of the pp force is determined by fitting the /3-decay 
half- lives of /3 emitters. This fitting procedure was systematically carried out in Ref. [8], where the strengths Xgt' 
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and kFqj, were considered to be smooth functions of the mass number A. The result found using a Nilsson potential 

as the deformed mean field was Xgt = ^-^ MeV and kFqj, = 0.58 /A^-'^ MeV. Nevertheless, this parametrization 
clearly depends on the model used for single particle wave functions and energies, as well as on the set of experimental 
data considered. Therefore, these coupling strengths can be used as a reference but cannot be safely extrapolated 
to different mean fields or different mass regions. As we shall see in the next section, the strengths from Ref. [8] 
reproduce well the data when using the WS potential, but one needs a somewhat smaller value of to reproduce 
the GT resonance with the HF mean field. 

We introduce the proton-neutron QRPA phonon operator for GT excitations in even-even nuclei 

^Ik = E + ^.T".",] , (2.1) 

where a+ (a) are quasiparticle creation (annihilation) operators, lok are the RPA excitation energies, and X':^^ ,Y^J^ 
the forward and backward amplitudes, respectively. The solution of the QRPA equations can be found solving first a 
dispersion relation [20], which is of fourth order in the excitation energies ook- 

In the intrinsic frame the GT transition amplitudes connecting the QRPA ground state |0) to one phonon states 
\u)k) satisfying 

r^, |o> = |o) = \ujk) (2.2) 

are given by 

{uJK\<JKt^\Q) = TM'^'' . (2.3) 

where 
and 



(2.5) 



where v's arc occupation amplitudes (?i^ = 1 — u^). 

Finally, the GT strength B{GT^) in the laboratory system for a transition /^/^^(O+O) — > IfKf(l^K) can be 
obtained as 



B{GT^) = 1^ \5Kfi (oJK \aot^\ O)^ + 25^,1 (ojk O)' 



(2.6) 



The 2i//3/3 decay is described in second order perturbation of the weak interaction as two successive Gamow- Teller 
transitions via virtual intermediate 1+ states. The half-life of the 2^(3(3 decay 



T'i%(0i-0i)= G'-^lMl'^^ 



(2.7) 



is given as a product of a phase space integral G"^" and the Gamow- Teller transition matrix element MqI^, which 
contains the nuclear structure effects. For a transition connecting initial and final ground states, it is given by 



0+ 1 1 OKt- \W7) {^7 1 ) i^K' 1 1 '^Kt- 1 1 0+) 



(2.8) 



where K = 0,±1 and mi,mf label the number of intermediate 1+ RPA states u>7 ^ ^7 reached from the initial 
|0^ > and final |0y > nuclear ground states, respectively. The overlap < i^7\'^7 > needed to take into account 
the non-orthogonality of the intermediate states reached from the initial ground state to those reached from the final 
ground state. It is given by [12] 



K^K'> = E 
d if 



Re,t, {BCSf\BGSi} , 



(2.9) 
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where £if.f label the quasiparticlc Triy pairs for the initial and final nucleus, respectively. The factor Refit includes the 
overlap of single particle wave functions of the initial and final nuclei [12] 




(2.10) 



The BCS overlap factor < BCSf\BCSi > is derived in Ref. [12]. An approximate expression is given by 



{BcsfiBcs,) « n n H< + n n (■ 





(2.11) 




with Nn the number of single particle states with the same values of parity and projection CI of the full angular 
momentum on the nuclear symmetry axis. 

We note that for the case of the spherical QRPA the derivation of the overlap factor of the intermediate nuclear 
states generated from the initial and final nuclei was outlined in Ref. [21]. For the case of the deformed QRPA the 
generalization of this derivation was presented in [12]. There, the importance of the BCS overlap factor, which is 
an integral part of the overlap factor of the two sets of the intermediate nucleus, on the evaluation of the double 
beta decay nuclear matrix elements was maintained. Then, this procedure of the derivation of the overlap factor of 
intermediate nuclear states was followed also in Ref. [15], however, with a significant approximation. The role of the 
BCS overlap factor was neglected. Of course, this can affect the final results significantly, especially in the cases when 
deformations of the initial and final nuclei are different. 



In this Section wo present results for the bulk properties of the nuclei under study based on the quasiparticle mean 
field description. We consider both WS potential and Skyrme HF approaches. 

In the case of HF, the first step is to study the energy surfaces as a function of deformation. For this purpose we 
perform constrained calculations [22], minimizing the HF energy under the constraint of keeping fixed the nuclear 
deformation. We can see in Fig. 1 the total HF energy plotted versus the quadrupole deformation parameter 



defined in terms of the microscopically calculated quadrupole moment Qp and charge root mean square radius rc- 

The results in Fig. 1 correspond to HF calculations with the force Sk3, which arc qualitatively similar to the results 
obtained with other Skyrme forces, such as SG2 and SLy4 [23] . We observe that the HF calculation predicts in some 
instances (Ge, Se, Zr, Cd, Sn, Te, Xe) the existence of two energy minima close in energy, giving rise to possible shape 
isomers in these nuclei. Solid lines represent the energy curves of the parent nuclei suffering the double beta decay, 
while dashed lines represent the energy curves of the corresponding daughter nuclei. 

We can see in Table 1 the experimental and the microscopically calculated charge root mean square radii rc with 
the forces Sk3 and SG2. Wc quote the two, oblate-prolate, results in those cases where the energies for the two shapes 
are very close. The values obtained for the charge radii are in good agreement with the experimental values from Ref. 
[24]. They are also in good agreement with the results obtained from relativistic mean field calculations [25]. 

In Table 2 we can see the theoretical and experimental quadrupole deformation. Experimental values have been 
extracted from the measured quadrupole moments from two different methods. In the first one the quadrupole 
deformation is obtained from Eq.(3.1), using the empirical intrinsic moments derived from the laboratory moments 
of Ref. [26] assuming a well defined deformation. In the second case the quadrupole deformations are taken from Ref. 
[27], where they were derived from experimental values of B{E2) strengths. In this case the sign cannot be extracted. 

Our theoretical values have been derived microscopically from the forces Sk3 and SG2, using the intrinsic quadrupole 
moments obtained as ground state expectations of the Q20 operator and the microscopic charge radii quoted in Table 
1. We show the results obtained for the equilibrium shapes using the Skyrme forces Sk3 and SG2. In those cases 
where a second minimum appears at close energy, we also show within square brackets the corresponding deformation. 
We compare our results with the results from relativistic mean field calculations of Ref. [25] and with results from 
systematic calculations [28] based on macroscopic-microscopic models (finite range droplet macroscopic model and 
folded Yukawa single particle microscopic model). The agreement between all the theoretical calculations is very 
remarkable and in general they are within the range of experimental values determined from Refs. [26] and [27]. The 
experimental values from Refs. [26] and [27] are represented in Fig.l by the endpoints of the black boxes. 



III. GROUND STATE PROPERTIES 




(3.1) 
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There is still another important experimental information relevant to double beta decay that we wish to explore here. 
This is the Qpp energy of the decay. The energy released in a double beta process in a transition from ground-state 
to ground-state is given by 

Q/3/3 = [Mparent - Mdaughter - 2me] = [2(m„ - TUp - TUe) + BE{Z, N) - BE{Z + 2,N-2)], (3.2) 

in terms of the nuclear masses M's, or similarly, in terms of the binding energies BE^s of parent (Z, N) and daughter 
{Z + 2,N -2) nuclei. 

We can see in Table 3 the experimental values of Qp/s- They are compared with the values calculated with the force 
Sk3. The agreement with experiment is reasonable taking into account that we are dealing with differences of energies 
ranging from 400 MeV in A=48 systems to 1200 MeV in A=150 and a tiny percent error in the theoretical masses 
may induce a relatively large error in (5/3/3 values. From Table 3 we can see that to match the experimental energies, 
we need to increase slightly the energy difference in A=48, 76, 82, 116, 128, 130, 136 and to decrease it in A=96,100,150. 

In this work we have used Sk3 as a representative of the Skyrme forces without any attempt to optimize the 
agreement with the experimental Q/3/3 values. Even if Sk3 force may not be the best Skyrme force to accurately 
predict the right (5/3/3 values, it is very interesting to see that experimental and theoretical values are in many cases 
quite close. For the cases examined here it appears that the (5/3/3 values are fairly well reproduced for A==82,100,150 
and that deviations from experiment in the worst case is 2.25 MeV in A=96, where (5/3/3(exp) = 3.35 MeV and 
(5/3/3(sk3) = 5.59 MeV. In the future it will be interesting to test other Skyrme forces which may give better fits to the 
(5,3ij values. A case of particular concern for double beta decay calculations is that of A=128 for which the Sk3 force 
gives Qfj/j = —0.10 MeV. Clearly, in this case further work is needed to make selfconsistent calculations that give the 
right (5/3/3 and GT strengths. 

We also notice that, as illustrated in Fig. 1, by changing slightly the deformation value in the vicinity of the equilib- 
rium deformation, one can change the binding energy correspondingly and get into agreement with the experimental 
Qpfj value. This procedure could be justified from the point of view that, in principle, in the HF method one could 
consider several collective degrees of freedom and that the absolute minimum in the multidimensional landscape could 
correspond to a slightly different /3 value. We can see in Fig. 1 the experimental value of 5 = BEi — BEf as the 
distance between the two horizontal lines plotted in each panel. The solid horizontal line refers to the energy of the 
parent while the dashed horizontal one refers to the daughter binding energy. One of them is always a reference and 
signals the energy to keep fixed (parent for A=96, 100,150 and daughter in A=48,76,82,116,128,130,136). The other 
line indicates the binding energy needed to reproduce the experimental S = BEi ~ BEf. Therefore, the cuts of this 
horizontal line with the corresponding energy curve indicates the deformations where this condition is satisfied. 

Fig. 2 shows the results obtained with HF-Sk3 for the GT strength distributions in ^^*Tc and ^■^^Xe, which are 
among the cases where the calculated (5/3/3 is worse and the change in /3 needed to fit (5/3/3 is larger. We show the 
results obtained with the equilibrium deformations (3 = —0.088 in ^^*Te and = 0.001 in ^^^Xe as well as with the 
deformations that fit the (5/3^3 values, f3 = —0.005 in -'^^^Tc and /? = 0.102 in ^"^^Xe. We can see that the strength 
distributions obtained with both deformations are similar except for a small displacement in energies. In all the other 
cases the effect is even smaller and the strength distributions obtained with the equilibrium deformation or with the 
slightly changed deformation are practically unchanged. In the next section we show GT strengths obtained at the 
HF minimum. 

In the case of the Woods-Saxon potential, where the deformation is an input parameter, we take the values from 
both Refs. [26] and [27]. Since for each nucleus, the two references give two different values of the /3— parameter, we 
show in the next section the GT distributions obtained with the two values to take into account this uncertainty. 

IV. GAMOW-TELLER STRENGTH DISTRIBUTIONS 

In this Section we show and discuss the Gamow- Teller strength distributions obtained from different choices of the 
deformed mean fields and residual interactions. 

We notice that the relevant strength distributions for the double beta decay are the B{GT~) distribution of the 
parent nuclei and the B{GT^) distribution of daughter nuclei. 

As a general rule, the following figures showing the GT strength distributions are plotted versus the excitation 
energy of daughter nucleus. The distributions of the GT strength have been folded with Gaussian hmctions of 1 
MeV width to facilitate the comparison among the various calculations, so that the original discrete spectrum is 
transformed into a continuous profile. These distributions are given in units of S'^/47r and one should keep in mind 
that a quenching of the gA factor, typically gA.cS = (0.7 — 0.8) .gA,frcc, which appears squared in the GT strength, is 
expected on the basis of the observed quenching in charge exchange reactions. 
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In the case of the B{GT^) distributions, we first observe the different scale, which is about one order of magnitude 
smaller than the B{GT^) scale. This is a consequence of the Pauli blocking. In the nuclei considered here the number 
of neutrons N is much larger than the number of protons Z . The difference between total B{GT~) and B{GT~^) 
strengths (Ikeda sum rule, which is fulfilled in our calculations), is then a large number given by 3{N — Z) and 
practically determined by the magnitude of the B{GT~) strength. 

We start in Fig. 3 with a discussion of the dependence of the GT strength distributions on the coupling strength 
of the particle- hole residual interaction Xgt ^'^^ ^ fixed value of the particle-particle coupling constant Kqj, = 0. The 
results correspond to HF with the force Sk3 in the A=150 case. We can see on the left panel the B{GT~) strength 
distribution of the parent nucleus ^"''"Nd and on the right panel the B{GT~^) strength distribution of the daughter 
nucleus ^^°Sm. The pairing gap parameters are given in Table 1 and the deformations are given in Table 2. We can see 
in Fig. 3 how the most important effect of Xcrr '^^ B{GT~) strength distribution is a shift of the strength toward 
higher excitation energies. This displacement of the GT strength is accompanied by a reduction of the strength. This 
reduction can be more clearly appreciated on the B{GT^) strength distribution because the scale in this case is about 
two orders of magnitude smaller, as it should be to fulfill the Ikeda sum rule YXB{GT-)-B{GT+)\ = 3{N-Z) = 90, 78 
for Nd and Sm, respectively. Therefore, the coupling constant Xgt pl^-ys an important role to reproduce the position 
of the GT~ resonance. On the other hand, the sensitivity of the GT strength distribution on the particle-particle 
coupling constant Kq'j, is not so important as it can be seen from Fig. 4, where we can see the GT strength distributions 
for a fixed value of X(jt = 0.156 MeV [8] and for several values of kJ^'!'^, on the example of HF with the force Sk3 in 
the A=150 case. As we can see, the position of the resonance does not change appreciably. Therefore, other methods, 
such as fitting the half-lives of unstable nuclei in the same mass region, have to be used to get phenomenologically 
their values. 

In the next set of figures (Figs. 5-14) we show, for each couple of double beta decay partners, the results obtained for 
the B{GT~) strength distributions of the parent nuclei on the top panels and for the B{GT^) strength distributions 
of the daughter nuclei on the bottom panels. Also shown are the experimental data whenever they are available. In 
each figure the left panels correspond to HF+BCS+QRPA calculations with the force Sk3 and the right panels to 
WS-I-BCS-I-QRPA calculations. In the case of HF we use the equilibrium deformations. We show with dashed lines the 
2qp results for HF-I-BCS calculations where the residual interaction is not considered. This serves as a reference and 
can be used to see the necessity of the residual force to get into agreement with experiment. Solid lines are the results 
obtained with x^r = 0-1 MeV and k^^, = &/A MeV, which produce a good fit to all the measured GT resonances of 
double beta emitters, as well as to the two-neutrino double beta decay matrix elements, as we shall see in the next 
section. This small value of the coupling constant needed to reproduce the experimental GT resonances within 
a selfconsistent approach with Skyrme forces is in agreement with the same observation reported in Refs. [11,13] and 
reflects the fact that one needs less residual interaction when using realistic effective density-dependent forces that 
when using phenomenological potentials to generate the single-particle energies and wave functions. 

In the case of calculations with the WS potentials shown on the right hand panels, we show results for the two 
different experimental deformations as obtained from Refs. [26] (solid lines) and [27] (dashed lines), which are given in 
Table 2. The calculations are done for a fixed value of the x^r ^^'^ '^gt constants as obtained from the parametrization 
in Ref. [8]. 

Some general common features to all figures can be established first. Concerning the HF calculations, the value 
of Xgt given by the parametrization of Ref. [8] is an overestimation when dealing with selfconsistent Skyrme HF 
calculations. Actually, a small value of Xgt = 0.1 MeV is already able to reproduce the experimental position of the 
GT resonance. This is a consequence of the structure of the two-body density-dependent Skyrme force that contains 
terms like spin exchange operators leading to a spin-spin interaction in the selfconsistent mean field, which is absent 
in the WS potential. The agreement with the experimental energy of the GT resonance is in this case very good as it 
can be seen in the cases A=76, 82, 100, 116, 128, 130, where this information is available. Indeed, the experimental giant 
GT resonances shown in these figures represent the centroids of broad bumps. The resonance in ^*Ca reported at 10 
MeV in Ref. [29] and used in the fitting procedure of Ref. [8] is also in good agreement with our results. 

With respect to the calculations performed with the WS potential, we can see that larger deformations produce 
peaks in the GT distributions displaced to higher energies. This is a consequence of the larger separation of the single 
particle energies when the deformation increases. Thus, since the deformation derived from Ref. [27] is larger than 
that of Ref. [26], solid lines appear in general on the left of dashed lines. 

It is also remarkable the good agreement with experiment obtained in this case. This agreement is mainly determined 
by the fixed value of Xgt from Ref. [8], which is still valid when describing the mean field with a WS potential. One 
should keep in mind that the parametrization of Ref. [8] was obtained using a Nilsson potential. 

In Table 4 wc; compare; the total GT strength measured and calculated with both HF and WS. When a standard 
quenching factor of 0.6 is included in the theoretical results fair agreement is found between theory and experiment. 
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V. TWO-NEUTRINO DOUBLE BETA DECAY MATRIX ELEMENTS 



In this section we analyze the effects of deformation, as well as the effect of the mean field and residual interactions 
on the 2u l3 l3-dccay. First, wc discuss the sensitivity to deformation of the nuclear structure contribution Mq^ to 
the 2vP(3 half-lives. This is done in Fig. 15, where we show the matrix elements Mq^^ as a function of both parent 
and daughter deformations. The figure corresponds to the decay ^^Zr — > ^^Mo calculated within a Woods-Saxon 
scheme with residual interactions from Ref. [8] . We have changed freely the deformations of both parent and daughter 
nuclei without any constraint from experiment. In this way we can study qualitatively the effect of deformation. The 
experimental values for Mq^^ shown in Fig. 15, as well in the next figures have been extracted from the adopted 
experimental half-lives T^j2 given in Ref. [34]. From the experimental half-lives and the corresponding kinematical 
factors G^'^ , we extract two experimental nuclear matrix elements from Eq. (2.7) by assuming values for the axial 
coupling constant qa = 1-25 or = 1. These two values are plotted in Figs. 15-18 as horizontal lines. 

From Fig. 15 we can see that the matrix elements have maximum values for equal deformations of both parent 
and daughter and these values decrease rapidly when the difference between the deformations of parent and daughter 
increases. In particular, we observe that the Mq^^ value obtained within a spherical picture (/3paront = /^daughter = 0) 
is about the upper limit and only comparable with values obtained with same deformations for parent and daughter 
in the deformed picture. As soon as the deformations of parent and daughter change, we get a reduction in the Mq^ 
matrix elements that cannot be obtained from a spherical description. The mechanism of this reduction due to the 
different deformations was studied in Ref. [12], where it was found that the overlap factor in Eqs.(2.8-2.11) is at the 
origin of the suppression. We can see in Fig. 15 that the experimental values of Mq^^, plotted as thick segments in 
each curve, are compatible with particular values of parent and daughter deformations. 

We show in Fig. 16 the difference between parent and daughter nuclear quadrupole deformations for the double-beta 
emitters. The dots correspond to the results obtained from selfconsistent HF calculations with the Skyrme force Sk3, 
while the extreme values on the vertical segments indicate the maximum and minimum differences compatible with 
the experimental values in Table 2. These are also the extreme values used in WS calculations in Fig. 17. 

In the next two figures (Figs. 17-18) we show the matrix elements Mq^ for the same double beta emitters studied 
in the last section as a function of the particle-particle strength Kq^. The experimental values extracted by assuming 
gA = 1.25 or gA = 1 are shown by the lower and upper horizontal lines, respectively. 

It is well known [4] that the pp interaction introduces a different mechanism of suppression of the Mq'^ matrix 
elements, which is also interesting to study in our case. In this way we can compare the effect of the pp force with 
the effect due to deformation. We first discuss in Fig. 17 the effects of deformation in the WS case by taking the 
available experimental quadrupole deformations and then we discuss the HF case by considering the selfconsistent 
deformations. 

In the WS case (Fig. 17) we use the same potential parameters, gaps and ph residual interaction as those used in 
the single beta calculations in the previous section. For deformations we take all the experimental possibilities for 
parent and daughter given in Table 2 and cross them to calculate Mq:^. Then we show in the figure the upper and 
lower results obtained as a function of n^lp and we draw a shadow region between them. We also show for comparison 
the results obtained in the spherical case (dashed lines). 

The first thing to notice is the already mentioned reduction of Mq^ as the magnitude of Kq^p increases, which takes 
place for both spherical and deformed cases, although the effect of the pp force is larger in the spherical case. The 
spherical curves decay faster than the deformed ones with Kq^. This means that the deformed results are more stable 
(more insensitive) to the particular strength of the pp interaction. 

Another interesting feature to mention is that, as expected from the analysis in Fig. 15, deformation introduces in 
most cases a reduction factor with respect to the spherical result. Only when the deformations of parent and daughter 
are very similar, the results obtained in the deformed case can be larger than the spherical ones. This is for instance 
the case of A = 82 (see Table 2), where the experimental quadrupole deformations [27] are /3 = 0.1944 in the case of 
the parent nucleus ®^Se and /3 — 0.2022 in the case of the daughter nucleus ®^Kr. 

In Fig. 18 we show the results corresponding to HF calculations with the force Sk3. In this case the quadrupole 
deformations are obtained selfconsistently (see Table 2) and are the same we used in the previous section to calculate 
the GT strength distributions in the single /?— decays. The coupling strength of the ph residual interaction has been 
taken x = 0-1 MeV as in the previous section. Contrary to the case of single /?— decay, where the position of the 
GT resonance is determined by the strength of the ph force and almost independent on the k^qj, force, we can see in 
Fig. 18 the sensitivity of the 2u/30 decay to the pp force. We find that the experimental values of Mq^ are roughly 
reproduced with values of k^j, = 6/A MeV. This is the reason why we used these values also in the HF calculations 
of the GT distributions in the previous section. 

To illustrate even further the effect of deformation on the M^^^ matrix elements, we show in Fig. 19 the 
HF-I-BCS-I-QRPA results with the Skyrme force SG2 for the decay ^^Zr — > ^^Mo corresponding to the decay of a 
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prolate parent to a prolate daughter or to the decay of a spherical parent into a spherical or a prolate daughter. The 
actual deformations used in these calculations are the consistently obtained with the force SG2 as can be seen in 
Table 2, namely (3 = 0.016, 0.147 for the parent ^^Zr and /3 = -0.006, 0.119 for the daughter ^^Mo. We can see that 
the transition from the spherical shape to the prolate shape reduces considerably the matrix elements as compared 
to the spherical/spherical or to the prolate/prolate cases, which are comparable. This reduction due to the different 
deformations makes the results compatible with the experimental values. 

VI. CONCLUDING REMARKS 

Using a deformed QRPA formalism, which includes ph and pp separable residual interactions, we have studied 
the GT strength distributions for the two decay branches and /?+ in double beta decay processes, as well as the 
two-neutrino double beta decay matrix elements. In the same manner in which two-neutrino double beta decay is 
used to calibrate the nuclear matrix elements for neutrinoless double beta decay, the single beta decay branches of 
parent and daughter are used to test the matrix elements for two-neutrino double beta decay. 

Two different methods, a deformed Skyrme HF approach and a phenomenological deformed WS potential, are used 
to construct the quasiparticle mean field, which includes pairing correlations in BCS approximation. 

In the case of HF the deformation is determined selfconsistently and we are able to reproduce the experimental 
charge radii and quadrupole moments. In the case of WS the input deformation is taken from experiment and we 
use two values for each nucleus, one corresponding to the measured quadrupole moment of the first 2"*" state and the 
other extracted from the measured B{E2) values. The latter can be considered as an upper limit of the /3— value. 
More experimental work would be needed to improve and complete the experimental determination of the quadrupole 
moments based on the first of these methods. 

Starting from these quasiparticle basis we perform a pnQRPA calculation with separable forces to obtain the energy 
distributions of the GT~ strength in the parent nucleus and the GT+ strength in the daughter nucleus and from them 
the 2i^/3/9— decay matrix elements. It is well known from previous studies that the ph force allows a reasonable fit of 
the GT resonance, and that the pp force also affects the B{GT~) and B{GT^) distributions. This knowledge has been 
applied in our paper to test existing parametrizations of the ph and pp residual forces and to see how good agreement 
can be obtained between experiment and theory with HF and WS methods. To our knowledge, we have considered 
for the first time simultaneously all the possible two-beta emitters and their corresponding daughters comparing their 
B{GT) to experiment. We find that we need different strengths of the ph force to reproduce the position of the GT 
resonance, depending on the HF or WS basis. In the first case, a small value of Xgt = 0-1 MeV reproduces all the 
measured GT resonances. In the second case the parametrization obtained in Ref. [8] {xgt ~ ^-^ /^°'^ MeV), using 
a Nilsson potential, is still valid when using a WS potential. The fact that Xgt smaller in HF than in WS can be 
understood as arising from the fact that the HF mean field already contains the average effect of spin-spin interaction 
terms. In both cases we reproduce reasonably well not only the position of the resonances but also the total GT 
strength. It should also be mentioned that the GT~ strength of the parent nucleus and the GT+ strength of the 
daughter are located at different energies, a feature that is relevant for double beta decay because it introduces a 
reduction of the double beta decay probabilities. It would be very useful to improve and complete the experimental 
information on GT strength distributions by (p,n) and (n,p) charge exchange reactions on nuclei participating in 
double beta decay. 

We have also explored the theoretical Q/3/3 values obtained with the HF method. We find that with the Sk3 force 
used here the agreement with the experimental Qi3p is not perfect and that it will be worth to look for a Skyrme 
force that may optimize agreement with experiment on both GT strengths and Q/3/3 values. However, taking into 
account that there is no fitting parameter at all, the agreement between theory and experiment is good. So far, no 
other approach used to study double beta decay has been able to obtain Qpp values so close to experiment. This is 
an important advantage of the HF method that deserves further exploration. 

The effect of deformation on the 2z//3/?— decay matrix elements has been studied first by considering the deformations 
of both parent and daughter as free parameters. It is found that the matrix elements are suppressed with respect to 
the spherical case. More precisely, it is found a sizable reduction effect that scales with the deformation difference 
between parent and daughter. This suppression mechanism, which is ignored in spherical treatments, may play an 
important role in approaching the theoretical estimates to experiment. 

In the case of the WS potential, we have studied the 2z//3/3— decay matrix elements by considering the maximum 
and minimum differences between the experimental deformations of parent and daughter. This procedure generates 
a region of decaying rates that would be narrowed from an improved experimental determination of the quadrupole 
deformations. 
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In the case of HF calculations we find that using the selfconsistcnt deformations obtained from the minimization 
of the energy and residual interactions with coupling strengths given by Xgt = ^-^ MeV and k^j, = 6/A MeV, we 
are able to reproduce simultaneously the available experimental information on the GT strength distributions of the 
single /3— branches and the 2i//3/3— decay matrix elements. 
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TABLE I. Pairing gaps for protons and neutrons Ap, A„ (MeV) and charge radii Vc (fm). 



Nucleus 




A„ 


exp Tc [24] 


rc Sk3 


Tc SG2 


rc [25] 


*«Ca 


2.18 


1.68 


3.4736(8) 


3.586 


3.549 


3.471 


48rj,j 


1.90 


1.56 


3.592 


3.628 


3.597 


3.571 


^«Ge 


1.56 


1.54 


4.127(8) 


4.130 


4.083 


4.057 


^«Se 


1.75 


1.71 


4.152(9) 


4.170-4.180 


4.113-4.143 


4.119 


«2Se 


1.41 


1.54 


4.122(8) 


4.204 


4.159 


4.131 


^^Kr 


1.72 


1.64 


4.1921(11) 


4.196 


4.196 


4.173 




1.53 


0.84 


4.3508(12) 


4.433-4.443 


4.342-4.389 


4.376 




1.53 


1.03 


4.377(10) 


4.448-4.457 


4.369-4.388 


4.381 


i°°Mo 


1.60 


1.36 


4.447(10) 


4.516 


4.439-4.466 


4.448 


^°°Ru 


1.55 


1.30 


4.453 


4.516 


4.457 


4.449 


"''Cd 


1.47 


1.37 


4.625 


4.703-4.715 


4.653 


4.643 


"«Sn 


1.77 


1.20 


4.625 


4.709-4.753 


4.702 


4.609 


128rpe 


1.13 


1.28 


4.735 


4.803-4.805 


4.746 


4.732 


i2*Xe 


1.32 


1.27 


4.776 


4.836-4.839 


4.782-4.786 


4.778 


iso^e 


1.06 


1.18 


4.742 


4.812-4.816 


4.750 


4.739 


i30Xe 


1.31 


1.25 


4.783 


4.845-4.846 


4.796-4.801 


4.784 


i36xe 


0.98 


1.44 


4.799 


4.878 


4.815 


4.804 


"«Ba 


1.27 


1.03 


4.833 


4.902 


4.847 


4.837 




1.23 


1.05 


5.047 


5.114 


5.055 


5.046 




1.44 


1.19 


5.047 


5.108 


5.046 


5.047 
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TABLE II. Theoretical and experimental /3 values, see text. 



Exp Theory 



NucloTis Rcf. [26] Rcf. [27] this work (Sk3) this work (SG2) Rcf. [25] Rcf. [28] 





0.000 


0.101(17) 


-0.002 


-0.001 


0.000 


0.000 


48rj,j 


+0.17(10) 


0.269(7) 


-0.002 


-0.003 


-0.009 


0.000 


^«Ge 


+0.095(30) 


0.2623(39) 


0.161 


0.157 


0.157 


0.143 


^«Se 


+0.163(33) 


0.3090(37) 


-0.181 [+0.157] 


-0.191 [+0.049] 


-0.244 


-0.241 


^^Se 


+0.104(32) 


0.1944(26) 


0.126 


0.150 


0.133 


0.154 


»2Kr 




0.2022(45) 


0.106 


0.103 


0.119 


0.071 






0.081(16) 


0.207 [-0.167] 


0.016 [+0.147] 


0.223 


0.217 


9«Mo 


+0.068(27) 


0.1720(16) 


0.147 [-0.164] 


-0.006 [+0.119] 


0.167 


0.080 


i°°Mo 


+0.139(30) 


0.2309(22) 


0.236 


0.167 [-0.191] 


0.253 


0.244 




+0.136(22) 


0.2172(22) 


0.175 


0.157 


0.194 


0.161 




+0.113(11) 


0.1907(34) 


0.206 [-0.207] 


0.209 


-0.258 


-0.241 


"«Sn 


+0.043(10) 


0.1118(16) 


0.264 [-0.134] 


0.251 [-0.034] 


0.003 


0.000 


128 rpg 


+0.011(10) 


0.1363(11) 


-0.088 [+0.102] 


0.094 [-0.091] 


-0.002 


0.000 


i28Xe 




0.1837(49) 


0.148 [-0.122] 


0.150 [-0.133] 


0.160 


0.143 


130 rpg 


+0.035(23) 


0.1184(14) 


-0.076 [+0.051] 


-0.039 [+0.066] 


0.032 


0.000 


i30Xe 




0.169(6) 


0.108 [-0.098] 


0.161 [-0.132] 


0.128 


-0.113 






0.086(19) 


0.001 


0.016 


-0.001 


0.000 


"6Ba 




0.1242(8) 


0.009 


0.070 


-0.002 


0.000 




+0.367(86) 


0.2848(21) 


0.266 


0.271 


0.221 


0.243 


i5°Sm 


+0.230(30) 


0.1931(22) 


0.207 


0.203 


0.176 


0.206 
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TABLE III. Experimental and theoretical Qfjp (MeV) values obtained with the Skyrme force Sk3. 



Double beta transition (Q,3,3)exp (Q/3/;)sk3 





4.272 


2.95 




2.039 


1.36 


82 se ^ ^^Kr 


2.995 


2.58 




3.350 


5.59 




3.034 


3.57 


ii6cd ^ "''Sn 


2.805 


1.88 


128Tg ^ 128xe 


0.867 


-0.10 


i30Te ^ "°Xe 


2.529 


1.20 


136xc ^ g,^ 


2.468 


n.80 




:i.:i()r 


:5.59 



TABLE IV. Experimental and calculated (HF-Sk3 and WS) summed GT strength. A standard quenching factor 0.6 has 
been included in the theoretical results. 



HF-Sk3 WS 



X;B(GT+) **Ti 1.42 ±0.2 1.00 1.79 

''®Se 1.45 ±0.07 0.48 2.06 

Y,B{GT-) '^^Ge 19.89 21.78 22.65 

^^Se 21.91 25.34 26.09 

i°°Mo 26.69 29.14 29.93 

"®Cd 32.70 34.79 36.41 

i2«Te 40.08 43.22 43.44 

"°Te 45.90 46.85 46.66 
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FIG. 1. Binding energy [MeV] £is a function of the quadrupole deformation parameter /3 obtained from deformed Hartree-Fock 
calculations with the Skyrme force Sk3. Experimental /3 values from Refs. [26] and [27] are represented as the extreme values 
of the black boxes. The experimental difference of binding parent and daughter binding energies are given by the distance 
between the two horizontal lines (see text). 
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FIG. 2. HF-Sk3 B{GT-) strength distributions [ffA/47r] in ^^*Te and ^^"Xe calculated with the equilibrium deformation 
(/3equii) and with the deformation that fits the experimental Qgg values (/3fit)- 
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FIG. 3. HF-Sk3 Gamow-Teller strength distributions [5i/47r] in ^^^Nd and ^^°Sni for various values of the coupling strength 
xfr [MeV]. 
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FIG. 4. HF-Sk3 Gamow-Teller strength distributions \g\l^-K\ in ^^^Nd and ^^°Sni for various values of the coupling strength 
K^^j, [MeV]. 
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FIG. 5. Gamow-Teller B(GT-) and B{GT+) strength distributions [5i/47r] in "^Ca and "^Ti plotted as a function of the 
excitation energies of the corresponding daughter nuclei. Left panels show results from HF(Sk3) calculations without residual 
interaction (dashed lines) and with residual interactions with Xgt ~ ^-lO MeV, k^j, = 6/ A MeV (solid lines). Right panels 
show results using WS potentials with Xgt ^'^'^ '^gt from Ref. [8] and with two different values for the quadrupole deformation: 
/3i from Ref. [26] (solid line) and P2 from Ref. [27] (dashed line). Experimental data are from Ref. [30]. Notice that no quenching 
factor has been included in the calculations. 
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FIG. 8. Same as in Fig. 5 for ^«Zr and ^^Mo. 
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FIG. 9. Same as in Fig. 5 for ^°°Mo and ^°°Ru. Vertical lines in ^°°Mo are experimental data from [33]. 
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FIG. 10. Same as in Fig. 5 for ^^®Cd and ^^®Sn. Vertical lines in ^^®Cd are experimental data from [33]. 
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FIG. 11. Same as in Fig. 5 for ^^*Te and ^^®Xe. Vertical lines in ^^®Te are experimental data from [32]. 
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FIG. 12. Same as in Fig. 5 for ^^°Te and ^^°Xe. Vertical lines in ^^°Te are experimental data from [32]. 
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FIG. 13. Same as in Fig. 5 for ^^^Xe and "®Ba. 
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FIG. 14. Same as in Fig. 5 for ^^°Nd and ^^^Sm. 
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FIG. 15. 2!^/3/3— decay matrix elements of Zr as a function of both parent and daughter deformations. The two dashed 
horizontal lines correspond to experimental Mqx extracted from Ref. [34] using qa = 1.0 and qa = 1.25. The thick segments 
in each curve correspond to the experimental values of Mq^- 
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FIG. 16. Difference between parent and daughter quadrupole deformations in double-beta emitters. Dots are selfconsistent 
results from Skyrme Sk3 calculations. Vertical lines indicate the maximum and minimum experimental differences (see Table 
2), which are used in WS calculations. 
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FIG. 17. 2i^/3/3— decay matrix elements calculated with Woods-Saxon potentials as a function of the particle-particle inter- 
action strength k^qj,. Dashed lines correspond to the results assuming spherical nuclei. Solid lines correspond to the results 
obtained by using the maximum and minimum differences between the experimental deformations of parent and daughter (see 
Table 2). Horizontal lines are the experimental Mqx extracted from Ref. [34] using qa = 1.0 and qa = 1.25. 



31 



exp 

HF-QRPA 



0.4 

0.2 


0.4 



> 0.4 



0.2 


0.2 

0.2 



1 1 1 1 1 1 1 

48 „ 48_. 

Ca -> Ti - 


^ 1 1 1 1 1 

- 76^ 76„ 

Ge^ Se 


- 


\ - 




- 


— 1 1 1 1 1 1 1 — 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 


^^^^Se Kr _ 


" 96_ 96, , \ " 

/ r ^ \An X 

Z-jL 7 IVHJ X 






— 1 1 1 1 1 1 1 — 


1 1 1 1 1 1 1 


1 1 1 1 1 1 1 


1 1 1 1 1 1 1 

116„, 116„ 

Cd^ Sn - 






_ 10(L ^ 100„ \ _ 

Mo ^ Ru \ 

— 1 1 1 1 1 1 1 — 




— ^ 

\ \ 1 1 \ \ 1 


1 1 1 1 1 1 1 

128„ 128„ 

""^.^ Te ^ Xe - 


1 1 1 1 1 1 1 

130_, 130„ 

Te — > Xe - 




— \ — \ — \ — \ — \ — I^H — 


— \ — \ — \ — \ — 1 P 1 — 

136,^ 136„ 

Xe ^ Ba . 


150-^, 150_ 

-^^^ Nd^ Sm. 






\ 1 


i\ 1 



0.4 

0.2 


0.4 

0.2 


0.4 

0.2 


0.2 



0.2 



0.02 0.04 0.06 0.02 0.04 0.06 0.08 



K [MeV] 
pp 



K [MeV] 
pp 



FIG. 18. 2J//3/3— decay matrix elements obtained from Skyrme(Sk3) deformed Hartree-Fock calculations as a function of the 
particle-particle interaction strength /t^. Horizontal lines are the experimental Mq^ extracted from Ref. [34] using = 1-0 
and QA = 1.25. 
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FIG. 19. 2J//3/3— decay matrix elements of ^^Zr obtained from Skyrme(SG2) deformed Hartree-Fock calculations as a function 
of the particle-particle interaction strength k^j-. The dotted curve corresponds to calculations using the prolate shapes for 

parent and daughter Pp — 0.147, /3d = 0.119 (see Table 2). The dashed curve is for spherical shapes and the solid curve is 
for spherical parent and prolate daughter f3p = 0.016, /3d = 0.119 (see Table 2). Horizontal lines are the experimental Mqx 
extracted from Ref. [34] using qa = 1.0 and qa = 1.25. 
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